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Abstract: The magnetoresistance (MR) in polycrystalline colossal magnetoresistive 
compounds follows a behavior different from single crystals below the ferromagnetic 
transition temperature. This difference is usually attributed to spin polarized tunneling 
at the grain boundaries of the polycrystalline sample. Here we derive a theoretical 
expression for the contribution of spin polarized tunneling to the magnetoresistance in 
ferromagnetic systems under the mean field approximation. We apply this model to our 
experimental data on the half metallic ferromagnet Lao.7Sr .3Mn03, and find that the 
theoretical predictions agree quite well with the observed dependence of the spin 
polarized MR on the spontaneous magnetization. 
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I. Introduction 

The large amount of recent activity on colossal magnetoresistive (CMR) 
manganites has revived interest in the study of electrical transport in granular itinerant 
ferromagnets. Many of these materials (e.g., Lao.7Sr .3Mn03, Cr0 2 , Fe 3 4 ) in the 
polycrystalline form, are now known to exhibit large magnetoresistance at low fields 
below the ferromagnetic transition temperature (r c ). u ' 3 ' 4 . Under similar conditions, 
magnetoresistance in single crystals is either very low or totally absent. The origin of 
this non-intrinsic intergranular effect is interesting and its detailed understanding 
should be important from the technological point of view. 

Though several researchers have attributed this effect to intergranular spin 
polarized tunneling, 1 ' 5 ' 6 ' 7 a complete understanding of this phenomenon is lacking. 
Comparing the magnetoresistance of single crystalline and polycrystalline samples of 
Lao.7Sr .3Mn03, Hwang et al. 1 were the first to suggest that spin polarized tunneling at 
the polycrystalline grain boundaries might play a crucial role in determining the 
magnetotransport properties below the ferromagnetic T c . They suggested that the 
magnetoresistance in polycrystalline LaovSro^MnC^ originates from the following two 
sources: (i) an intrinsic part arising from Zener double exchange mechanism 8 between 
two neighboring manganese ions, and (ii) the intergranular spin polarized tunneling. 
The second component produces a sharp drop in resistance at low fields and is 
dominant at temperatures much below T c , whereas the first component is dominant 
close to T c . In an earlier paper 6 we have shown that the intergranular part in the 
magnetoresistance can be distinguished from the intrinsic part by considering the field 
dependence of the magnetoresistance. The contribution from the intergranular spin 
polarized tunneling gives rise to a sharp drop in resistance at low fields and low 
temperatures. At high fields, the resistance varies almost linearly with the field. The 
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high field behavior can be attributed to the intrinsic Zener double exchange mechanism 
and explains the observations in single crystals. Thus one can estimate the resistance 
drop due to spin polarized tunneling by finding the intercept obtained by extrapolating 
back the linear high field region. Using this scheme, we had earlier shown that the 
intergranular magnetoresistance due to spin polarized tunneling in polycrystalline 
colossal magnetoresistive manganites drops monotonicalfy with temperature, whereas 
the intrinsic part follows the behavior expected from Zener double exchange. 6 

In an attempt to understand the mechanism of magnetoresistance in granular 
ferromagnets, Hellman and Abeles 9 had proposed a model based on intergranular spin 
polarized tunneling. This model assumed that when an electron tunnels across a 
boundary between two grains having antiparallel magnetization, it will experience a 
potential barrier of the order of the exchange energy (E m ). So the tunneling probability 
between the two grains with antiparallel magnetization was assumed to be reduced by 

a factor e~ E,nlkT from that when they are parallel. The conceptual difficulty in 
Hellman and Abeles' model is the following. When the conduction band in the 
ferromagnet is partially polarized, the dominant factor guiding the tunneling probability 
is the respective up and down density of states (DOS) in the two grains. Inoue et al. lQ 
had pointed out the difficulty in taking an exponential factor in which the energy is of 
the order of exchange energy. Thus, the mechanism of spin polarized transport in 
granular ferromagnets needs to be investigated in greater detail. 

In this paper, we propose a model to understand the mechanism of spin 
polarized tunneling in granular ferromagnets. Section II describes the theoretical 
model, while the sample preparation and experimental procedures are given in Section 
III. In Section IV, we compare the predictions of the theoretical model with the 
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temperature and field dependence of spin polarized tunneling in Lao.7Sro.3Mn0 3 . We 
also discuss the applicability of this model to other granular ferromagnets. 



II. Theoretical model 

In a granular ferromagnet in zero field, the magnetization directions of the 
grains are randomly oriented due to the random orientation of their magnetocrystalline 
anisotropy axes. Following Inoue et al., 10 we first derive an expression for the 
difference in the spin polarized tunneling resistance between a configuration where the 
magnetization of two grains are at an angle 0, and that in which their magnetizations 
are made mutually parallel by the application of a magnetic field. Let nr and nj, be the 
respective densities of states of up and down spins at the Fermi level. We choose our 
z-axis to be parallel to the magnetization direction of grain 1 (see Fig. 1). The up and 



down spin eigenstates are defined as | S z ; T) = and | S z ; i) = 



0) 



respectively. 



We define the spin eigenstates in grain 2 as | S 9 ; T) and | S 9 ; -l). Using the Pauli 
rotation matrix R(0), we can easily find the relation between | S e ; T-l) and | S z ; T-l) 
from 

^cos(0/2) -sin(0 /2j\ 



I S 9 ; U> = R(0) | S z ; U>, where R(0) 
Hence, we obtain 



sin(0 / 2) cos(0 / 2) 



(1) 



I S 9 ; T> = cos(0/2) | S z ; T> + sin(0/2) | S z ; i), 
S e ;i) = -sin(0/2) | S z ; T> + cos(0/2) |s z ;i>. 



(2) 



4 



An electron can tunnel from grain 1 to grain 2 using one of the following channels: 



(a) 


|S Z ;T> -> 


1 S ; T> 


(b) 


|s z ;T>^ 


Is e ;i> 


(c) 


|S Z ;I>^ 


|S 9 ;T> 


(d) 


1 S z ; i) -> | 


S e ;i>. 



Since the Hamiltonian involved with the tunneling process is spin independent, the 
matrix elements corresponding to (a) and (d) should be proportional to cos 2 (0/2), and 
the matrix element corresponding to (b) and (c) should be proportional to sin 2 (0/2). 
The total transition probability of up (or down) spin for tunneling from grain 1 to grain 
2 should also depend on the initial and final DOS available to it in the two grains. 
Thus, the total transition probabilities for the four processes will be 

T a oc riT 2 cos 2 (0/2) 

T b °< nr nisin 2 (0/2) 

T c nr nisin 2 (0/2) 

T d oc m 2 cos 2 (0/2). 

The tunneling conductivity, o(0), between the two grains at angle involve the sum of 
all four processes. Thus, 

o(0) oc (nr 2 +ni 2 )cos 2 (0/2) + 2nr n|sin 2 (0/2) 
or, G(0) oc (l/2)(nr+ni) 2 [l+P 2 cos0], (3) 
where P = (nr-nj,)/(nt+nj,). When the magnetizations of grain 1 and grain 2 are made 
parallel to each other by applying a magnetic field, the conductivity becomes: 

0(0=0) oc (n T 2 W). (4) 
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In a real system with many grains and grain boundaries one has to average over cosG. 
Thus, the resistance change (AR spt ) arising due to spin polarized tunneling is 
proportional to [l/a(9)-l/a(0)]. 

In order to determine the temperature dependence of AR spt , we require the 
temperature dependence of nt and nj,. We use the ferromagnetic Kondo Hamiltonian 
initially proposed for this system by Furukawa 11 ' 12 ' 13 to calculate the up and down 
DOS. 

H = ~tH Cia Cja ~ Jn'LOi- Si 

(ij) i 

where t is the nearest neighbor hopping energy of the e g electrons and J H is the local 

S 

ferromagnetic Hund's rule coupling between the e g electron spin , and the ? 2g spin o; at 
the i-th site. In terms of a ± = o x ± z'o y , o z , S* = S x ± is y we can write, 

H = -t I cL cja -^L(otsi + at sT + 2 of S f ) 

(i.J) 2 i 

We resort to a mean field approximation in which the f 2g spin a is treated as a classical 
vector, i.e., the operators cr and (f are replaced by c-numbers: 

ar = A . ct = A . and Gi = B • 

We also write the e g spin operators in terms of annihilation and creation operators: 

st = Cft Cu , Si = Cu Cft and Si = \ (cj t Cft _ cjj, CfJ • 
As a result, we get the following expression for the mean field Hamiltonian: 

H MF = ~t L Cia Cja ~ ^ L [A Cft Cn + A* Ca Cft + B (c- t Cft ~ cj; CfJ] • 



6 



Going over to A:-space by Fourier transforming according to 



t ik.r 



t _ V t 

Cf C -Lc ka e 
k 



we get, 



h- = E [(- r„ - \ b y„)ci T c*t + (- r* + \ b y.JcLcu - y (AcI T cw + A*4i ct) 



Here = *Xe*" A where A is the nearest neighbor vector. We can rewrite this as 



h mf = l( 



* \Ckl) 



where 



h k = 



This 2x2 matrix is easily diagonalized to get the eigenvalues as 



kt = U-2y t ±j H a). 



i_ 2 — 2i 2_i_ 2 — 

where a - o x + O y + O z = 



2 ■ 2 



+ g . The corresponding eigenvectors are: 



- Ai) and \ Uk j- v. 



Uk — Xk 



where *J = ^ (l + f)and ^ = 1 (l-f). and. | T) = (>)a„d 1 1) = (J) . 



At zero temperature the ? 2g spin assumes a saturation value, so Bio = 1, and as 



a result the ground state 



) has purely up-spin character. The density of states 



corresponding to T=0 is shown schematically in Fig. 2(a). For T > 0, B/o < 1 and as a 
result the ground state and excited state have mixed character: 



Nt(E) = x\ N(E), andNi(E) = y\ N(E). 
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The up and down spin characters, xl and y\ , respectively, are both non-zero in this 

case, as shown in Fig. 2(b). Finally, for T>T C , B vanishes, resulting in x\ = y\. This 
leads to N(E) being an equal admixture of up and down spins 14 as shown schematically 
in Fig. 2(c). In our theory, the magnetization of the e g spin is: 

Therefore the total (spontaneous) magnetization is M = (B/2o) + B = aB, where a = 
(1 + l/2o). We can thus write xl and y\ in terms of the reduced magnetization, m = 
M(T)/M(T^0) as: 

xl = (1 + m)/2, and y\ = (1 - m)/2 (5) 
Within our theory, the ratio Nt(E)/Ni(E) is independent of energy. Thus, the up 
and down DOS at the Fermi level, nt and ni, which appear in equations (3) and (4) are 
proportional to xl and y\ respectively. Substituting equation (5) in the expression for 
spin polarized tunneling conductivity in equation (3) and (4), we obtain a relation 
between AR spt and the reduced spontaneous magnetization of the system: 

Atf spt oc[<l/ O (0)>- (l/o(0)] (6) 
Here the conductivity in zero field, l/o(0), is averaged over different values of cosG for 
various grains. While fitting AR spt vs. m data (see Section IV) we have used (cos(0)> as 
a fitting parameter. 

III. Experiment 

Bulk samples of La .7Sro. 3 Mn03 were prepared by a wet chemical route. 
Initially, a stoichiometric mixture of La 2 3 and Mn was dissolved in nitric acid and 
Sr(N0 3 )2 added to it. The cations were co -precipitated as carbonates by adding excess 
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ammonium carbonate to the nitrate solution. This carbonate precursor was then heated 
at 1200 C to obtain single phase Lao.7Sr .3Mn03. Electrical measurements were carried 
out on a rectangular bar shaped sample (2mm x 3mm x 15mm) using the conventional 
four probe technique. Magnetoresistance measurements were made using a home made 
superconducting magnet up to fields of 3 Tesla, in the temperature range 5 K to 300 
K. Magnetization versus field isotherms were measured on a Quantum Design SQUID 
magnetometer up to a field of 5.5 Tesla. 



IV. Results and Discussion 

Figure 3(a) shows the magnetoresistance, MR = [R(H)-R(0))/R(0], as a 
function of the applied field (H) from 5 K to 300 K. We observe a sharp drop in the 
MR at low fields and the magnitude of the drop decreases with increasing temperature. 
This sharp drop is absent in single crystals and epitaxial thin films 1 ' 2 ' 6 and is therefore 
associated with the grain boundaries in polycrystalline samples. At high field, the slope 
of this curve becomes almost linear. The slope increases with increasing temperature 
and is identical to the nature of the MR vs. H curves observed in single crystals. This 
part of the MR can be ascribed to the Zener double exchange and arises from the 
intrinsic magnetoresistance of the grains. 6 To find the total resistance drop associated 
with the intergranular part (A/? spt ) we extrapolate back the linear high field region to 
find its intercept at zero field. This is schematically shown in Fig 3(b) for the curve 
obtained at 5 K. Figure 3(c) shows the magnetization versus field (M-H) data for the 
same sample from 5 K to 340 K. The inset shows the spontaneous magnetization (M s ) 
as a function of temperature calculated by back extrapolating the linear region of the 
curve beyond technical saturation. 
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In Fig. 4 we have plotted Ai? spt as a function of m ( = M S (T)/M S (5K) ) along 
with the theoretical curve (Eq. 6). The proportionality constant in Eq. (6) and (cosG) 
are taken as fitting parameters. We get the best fit for (cosG) = -0.8788, which 
corresponds to 0-151.5 . The negative value of (cosG) is consistent with the fact that 
there are two interactions governing the direction of magnetization of a grain: the easy 
axis which is random and tends to randomize the direction of the magnetization, and 
the magneto static dipolar interaction between the grains which tend to favour 
antiparallel alignment between grains. 

So far, we have been concerned with the total MR arising from spin polarized 
tunneling, that is the total drop in resistance due to spin polarized tunneling after the 
technical saturation in magnetization. We now investigate the MR versus M at low 
fields. This is the regime in which domain rotation takes place, i.e. where (cos(0)) is 
changing. The MR versus (M/M s ) curves at fields less than 1 kOe are shown in Fig. 6. 
These curves fit well with a power law behavior: MR (M/M s ) n , where n<l. 
However, to fit these curves within our model, we need to compute (cosO) as a 
function of field, and this is beyond the scope of the present work. 

We have thus shown that the model of spin polarized tunneling proposed in this 
paper fits well with the data on the half metallic ferromagnet La .7Sro.3Mn03 (with the 
underlying t 2g spins forming a localized magnetic lattice and the e g spins forming a 
polarized conduction band). However, equations (3) and (4) should apply equally well 
to other granular ferromagnetic systems, provided the evolution of the up and down 
spin density of states with temperature is known. It might therefore be also interesting 
to investigate the nature of spin polarized tunneling in more conventional granular 
itinerant ferromagnets within the realms of this model. 
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V. Conclusion 

In summary, we have developed a mean field model for spin polarized 
tunneling in granular itinerant ferromagnets. The theoretical model provides a 
satisfactory fit for our data on the temperature dependence of the magnetoresistance 
due to spin polarized tunneling in the half metallic ferro magnet Lao.7Sro.3Mn0 3 . Further 
investigations on other systems should provide a better understanding of the transport 
mechanism in itinerant magnetic systems. 
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FIGURE CAPTIONS 



FIGURE 1. Two neighboring ferromagnetic grains with magnetization Mi and M 2 
making an angle with respect to each other. In a sufficiently large applied field, 
becomes zero, thus increasing the tunneling probability of an electron from grain 1 to 
grain 2. In our convention, the z-axis is along the magnetization direction of grain 1 . 

FIGURE 2. A schematic diagram of the up and down densities of states at various 
temperatures: (a) T=0, (b) 0<T<T c and (c) T>T C . E F shows the Fermi energy for the 
hole doped compound. At finite temperature, the spin eigenstates have a mixed 
character. 

FIGURE 3. (a) Magnetoresistance versus field isotherms for polycrystalline 
Lao.7Sr .3Mn0 3 from 5 K to 300 K. The sharp drop at low fields arises from spin 
polarized tunnelling at the polycrystalline grain boundaries, (b) The total resistance 
drop due to spin polarized tunnelling (A/? spt ) is estimated by back extrapolating the 
linear high field slope and finding its intercept at zero field, (c) The magnetization 
versus field isotherms at different temperatures. Inset shows the temperature variation 
of spontaneous magnetization (M s ) calculated by back extrapolating linear portion 
curve after technical saturation. At 340 K the spontaneous magnetization is calculated 
from the linear portion of the Arott plot (M/H versus M 2 ). 
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FIGURE 4. The total resistance drop due to spin polarized tunneling (Ai? spt ) at 
different temperatures as a function of the reduced magnetization m = M S (7)/M S (5K) 
(filled circles). The solid line is the theoretical best fit curve (Eq. (6)) to the 
experimental points. The best fit value of <cos(0)> is -0.8788. 

FIGURE 5. Magnetoresistance versus magnetization isotherms at low fields (up to 
1 kOe). The solid lines are the best fit to the relation MR (M/M s ) n . Note that this 
relation provides a good fit to the data over the entire temperature range. 
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